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The excited state dynamics of correlated electron and electron-phonon systems triggered 
by an oscillating electric-field pulse of large amplitude are theoretically investigated. A 
“negative-temperature” state and inversion of electron-electron and electron-phonon inter¬ 
actions are induced even by a symmetric monocycle pulse. This fact is numerically demon¬ 
strated, using the exact diagonalization method, in a band-insulator phase of one-dimensional 
three-quarter-filled strongly dimerized extended Peierls-Hubbard and Holstein models. When 
the total-energy increment is maximized as a function of the electric field amplitude, the occu¬ 
pancy of the bonding and antibonding orbitals is inverted to produce a negative-temperature 
state. Around this state, the dependences of time-averaged electron-electron and electron- 
phonon correlation functions on interaction parameters are opposite to those in the ground 
state. 


1. Introduction 

Photoinduced changes in electronic properties of correlated electron and electron-phonon 
systems have been investigated extensively so far . 1 - 2 Development in experimental techniques 
has allowed for progressively shortening optical pulses and consequently improving time 
resolutions. At the same time, the field amplitude that can be applied to materials has sig¬ 
nificantly increased for different excitation energies. This fact enhances the possibility of 
controlling the electronic properties. 

As for events that occur with an oscillating electric field of large amplitude, dynamical 
localization is known in the continuous-wave case . 3,4 When the frequency of the field is com¬ 
parable to or higher than that of electron transfers determined by the transfer integral, the slow 
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component of the electron dynamics is governed by the time average of the transfer integral 
with the Peierls-phase factor. 5 The time average for a sinusoidally oscillating field gives an 
effective transfer integral modulated by the zeroth-order Bessel function whose argument is 
proportional to the ratio of the amplitude to the frequency of the field. Thus, the effective 
transfer integral can vanish to cause dynamical localization or can even change from positive 
(negative) to negative (positive). 

Applying such a strong field to interacting electron systems, Tsuji et al. have shown that 
the sign inversion of the transfer integral and the consequent inversion of the band structure, a 
“negative-temperature” state, are regarded as the sign inversion of the interaction in the half- 
filled Hubbard model, using the dynamical mean-field theory, by which a superconducting 
state is possibly induced. 6 In addition, they have shown that such a repulsion-to-attraction 
conversion is realized not only by continuous waves but also by half-cycle pulses and asym¬ 
metric monocycle pulses. 7 In the limit of the vanishing pulse width, the irradiation causes 
a sudden change in the Peierls phase of the transfer integral, and the dynamical phase shift 
is equivalent to the crystal momentum shift. In this limit, however, a symmetric monocy¬ 
cle pulse causes nothing. Indeed, the repulsion-to-attraction conversion is hardly realized by 
symmetric monocycle pulses with various trial time profiles of finite pulse widths. 7 

The effective transfer integral modulated by the application of a continuous wave and 
that by the application of a pulse are not always the same. It depends on the system and 
how the system is excited. The quasi-two-dimensional metal complex Et 2 Me 2 Sb[Pd(dmit ) 2 ]2 
(dmit = l,3-dithiol-2-thione-4,5-dithiolate) is known to show a photoinduced transition from 
a charge-ordered-insulator phase to a Mott-insulator phase. 8 " 10 In this particular case, it has 
been demonstrated that continuous-wave lasers modulate the effective transfer integrals be¬ 
tween the dimers, and pulsed lasers modulate that within the dimer, both through the zeroth- 
order Bessel function. 11 On the basis of this fact, the dynamical localization is transiently 
expected for pulsed lasers if the field amplitude is sufficiently large. Indeed, optical freezing 
of charge motion has been found in an organic conductor, ar-(BEDT-TTF) 2 l3 (BEDT-TTF = 
bis[ethylenedithio]-tetrathiafulvalene). 12 Although its mechanism is not theoretically under¬ 
stood yet, it is regarded as due to collaboration with the dynamical localization and inter- 
molecular charge correlations. 

The realization of a negative-temperature state or modulation of transfer integrals and 
interaction parameters requires an oscillating electric field of large amplitude. In general, it 
is much larger than what is needed for some phase transition to be simply induced. In any 
case, the fact that transfer integrals are effectively modulated by pulsed lasers is encour- 
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aging because the field amplitude can be increased by shortening the pulse width. In the 
actual case of a'-(BEDT-TTF) 2 I 3 above, 12 the dynamical phase shift is known to be much 
smaller than n. On the basis of these findings, we study the effect of symmetric monocy¬ 
cle electric-field pulses whose time integrals are zero, i.e., whose dynamical phase shifts are 
zero. It was expected to be difficult to form a negative-temperature state or to invert inter¬ 
action parameters by such symmetric pulses. 7 The model systems we treat are similar to 
that used in Ref. 11, but they are simplified to use the exact diagonalization method: one¬ 
dimensional three-quarter-filled strongly dimerized extended Peierls-Hubbard and Holstein 
models. Their band-insulator phases consist of almost half-filled dimers alternating with al¬ 
most completely filled dimers. Charge transfers take place between these dimers, but the 
optical gap is determined by the electronic transition within the almost half-filled dimer, as in 
Et 2 Me2Sb[Pd(dmit) 2 ]2. 

In this paper, we show that, even by a symmetric monocycle electric-field pulse, a 
negative-temperature state is produced, and interactions are effectively inverted when we 
interpret the dependences of the time-averaged correlation functions on the interactions. 
Namely, the stronger on-site repulsion causes the larger space-time-averaged double occu¬ 
pancy, the stronger nearest-neighbor repulsion causes the larger averaged nearest-neighbor 
charge-density correlation, and so on. Such a negative-temperature state appears when the 
total-energy increment is maximized as a function of the electric field amplitude. 


2. Strongly Dimerized Models with Charge-Ordered Ground States 

We study the electron-phonon dynamics after strong photoexcitation in a one-dimensional 
three-quarter-filled strongly dimerized model with on-site and nearest-neighbor repulsions 
and with the Peierls and Holstein types of electron-phonon couplings, 

N/2 

H = ~ ^|(fi> ~ U j) ^i ( C 2j-1'Q- C 2j,cr + C 2j,Q- c 2j-l,cr) 


7=1 o- 

N/2 


X /mter Xi (' 


2 7 > 


c 2j+l,, 


+ C 


t 

2j+l,cr 



7=1 cr 


N N 

+ u ^ « ( ,T n U + V ^ n,n i+l 

i= 1 i= 1 


Z 


‘ u ± 

2S 


+ 


2 So) 




l) 


N / 3\ N 

+§ Y^ bi + b i) [ n i - 2 + 2 b 'i bi ’ 

r=i ' f i=i 


(1) 



J. Phys. Soc. Jpn. 


(a) 


i=l, 2, 



Transfer: t i ^inter 

= t 0 —n l 

Distance: Cl j a inter 


f 2 (liter 

= ? 0 — w 2 

^inter 



Fig. 1. (a) Schematic representation of transfer integrals and distances between neighboring sites, (b) Half- 

filled dimer whose bonding orbitals are occupied and completely filled dimer in ground state, (c) Half-filled 
dimer whose antibonding orbitals are occupied and completely filled dimer in negative-temperature state. 


where c] (j creates an electron with spin cr at site i, n i lT = cj^c^, and n, = X, r «,>• The transfer 
integral within the j-th dimer (i.e., between sites 2 / - 1 and 2 j) is modulated by the lattice 
displacement uf. t 0 - uj. The transfer integral between the dimers is denoted by f inter . They are 
shown schematically in Fig. 1(a). The parameter U represents the onsite repulsion strength, 
and V represents the nearest-neighbor repulsion strength. The lattice displacement Uj has a 
bare frequency co P and is coupled with the bond density between sites 2j - 1 and 2 j with 
strength S . It is noted that u ’s have a dimension of energy. The operator b] creates a phonon, 
which has a bare frequency <x> H and is coupled with the charge density relative to its average 
value at site i with strength g. 

The periodic boundary condition is imposed. For g=0, we use iV=16 sites. For g ± 0, we 
use iV=8 sites and limit the number of phonons up to a maximum of two on every site. The 
number of electrons is 24 for N=\6 and 12 for N=8 at three-quarter filling. The electron- 
quantum-phonon system is treated by the exact diagonalization method. The lattice system 
is treated classically. This is because the lattice vibration is caused by the intermolecular 
restoring force and consequently is slow compared with other dynamics. Meanwhile, the 
molecular vibration is caused by the intramolecular restoring force and its energy is often 
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comparable to those of the intermolecular charge transfers. Its quantum nature is important 
and essential for the polaron effect. Before the photoexcitation, the system is in the lowest- 
energy state so that the lattice displacements are determined by 


- S ^ (( c 21-1^2j,cr) + ( C 2j,a- C 2j-I,cr)) ■ 


( 2 ) 


They spontaneously alternate large and weak distortions, to which u\ and w 2 are assigned 
[Fig. 1(a)]. The electron-quantum-phonon ground state is in the band-insulator phase, which 
is schematically shown in Fig. 1(b) with q > t 2 > 0: it consists of strongly ( u\ < 0, large 
| Mi |) distorted dimers with the bonding orbitals occupied and weakly (m 2 < 0, small | w 2 I) 
distorted dimers with both the bonding and antibonding orbitals occupied. 

Photoexcitation is introduced through the Peierls phase 
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where r, 7 = r 7 - r, and r, being the location of the i-th site, which is substituted into Eq. (1). 
For symmetric monocycle electric-field pulses, we use the time-dependent vector potential, 
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where F is the field amplitude and the central frequency io is chosen to be nearly resonant 
with the optical gap. The time-dependent Schrodinger equation for the exact many-electron- 
quantum-phonon wave function is numerically solved by expanding the exponential evolution 
operator with a time slice dt- 0.02 to the 15th order and by checking the conservation of 
the norm. 1314 We use the leapfrog method to solve the classical equations for the lattice 
displacements, 

1 dru: 


Scoj, dt 2 


Fj({u}) , 


where the force Fj({u}) is obtained with the aid of the Hellmann-Feynman theorem: 
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( 6 ) 


For the bare transfer integrals, we use t 0 =l and t inter =0.1, which correspond to an extrin- 
sically and strongly dimerized system. The interaction parameters U, V, g, and S are varied. 
The Holstein coupling strength g is chosen to be so weak that the number of phonons is always 
much smaller than one on every site, ('P(t) | Id b; \ V F(/)) <sc 1. The bare phonon frequency co H , 
corresponding to a molecular vibration frequency, is set to be o) H = 0.2. The bare phonon fre¬ 
quency cl>p, corresponding to an intermolecular oscillation frequency, is set to be a> P = 0.05 if 
only S is nonzero among the couplings U, V, g, and S (i.e., in the Peierls model) and cl> p = 0 
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Fig. 2 . (Color online) Energy increment and Jo(eaF/Hcj) sin (eaF/tico) term for monocycle nearly resonant 
(co = 1.05 wq) excitation of one-electron two-site tight-binding model. 11 


with S = 0.1 otherwise (i.e., in the extended Hubbard and Holstein models with frozen lattice 
displacements). The intermolecular distances r (/ - are fixed for simplicity and so chosen that 
the Bessel functions Jo{erijF/fuo) modulating the transfers between sites i and j = i + 1 0=1, 
3 for the two nonequivalent intradimer transfers and i= 2, 4 for the interdimer transfers) are 
easily distinguishable: a\ = r l2 = a, a 2 = r 34 = 1.2a, and ai nter = r 2 3 = r 45 = 1.6a. For 
S = 0.1, q = to - Mi - 1.2 and t 2 = t 0 - u 2 - 1.0 [Fig. 1(a)]. The difference between t\ and t 2 
becomes smaller for larger U. 

As pointed out previously in Appendix of Ref. 11 for the two-site tight-binding model 
and for a resonantly excited strongly dimerized model, the total-energy increment AE lot is 
roughly proportional to J 0 (eaF/tuo) sin(eaF/Tito) although A£ lot is nonnegative. For refer¬ 
ence, we show again, in Fig. 2, the energy increment [Eq. (A-ll) in Ref. 11] and the term 
proportional to J 0 (eaF/fia>) sin(eaF/fitL>) [Eq. (A-13) in Ref. 11] for the monocycle nearly 
resonant (co = 1.05o» 0 with the bonding-to-antibonding transition energy frto 0 = It) excita¬ 
tion of the one-electron two-site tight-binding model with transfer integral t (upper panel of 
Fig. A-3 in Ref. 11). The maximum of J 0 (eaF/hto) sin (eaF/fLco) appears around eaF/hco= 1, 
while the maximum of the energy increment appears at eaF/tico slightly smaller than 1. 

3. States after Symmetric Monocycle Pulse Excitation 

As a typical case, the time evolution of the charge density ( v P(f) | /7, | V F(/)) and that of 
the double occupancy ('Fit) \ \ T(/)) are shown in Figs. 3(a) and 3(b), respectively, 

for U = 0.15, V = 0, g = 0, and 0 < t < 50 T. The charge density is, initially and a short 
time after the pulse excitation t > 5T, almost 1, 1, 2, 2 for sites i= 1, 2, 3, 4, respectively 
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Fig. 3. (Color online) (a) Charge density ('¥(t) | n, | 'P(f)) and (b) double occupancy ( v P(f) | n/^n/.x I v P(t)>, 
at four nonequivalent sites during and after monocycle pulse excitation with eaF/tiij =1.0 and <x> = 2.4. The 
model parameters are U = 0.15, V - 0, g = 0, S =0.1, and top = 0. The time averages below are taken for the 
time region indicated by t w . 

[see also Fig. 1(a), denoted formally by the (1,1,2,2) configuration hereafter]. Because the 
present system is in the band-insulator phase, the time evolution of the double occupancy 
is quite similar to that of the charge density. The time-averaged double occupancy below is 
calculated by 



(7) 


■W Jt s 


with t s = 5 T,t w = 45 T, and T = 2n/a> being the period of the oscillating electric field. The 
time domain of integration is indicated by the two-headed arrow. The other time-averaged 
quantities are calculated likewise. 

3.1 Effect of on-site repulsion U 

Time-averaged quantities for V = 0 and g - 0 with different values of U and fixed a>=2A 
are shown in Fig. 4 in the range of 0 < eaF/ha> < 2.5. The total-energy increment AE tot is 
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Fig. 4. (Color online) (a) Total-energy increment AF tm , (b) time-averaged double occupancy averaged over 
strongly distorted dimer sites (((ni,|Hi,|)) + ((« 2 ,t w 2 ,j .)))/2 and that averaged over weakly distorted dimer sites 
(((»3,T H 3,l)) + (( n 4,T w 4,i)))/2- (c) average double occupancy 4 Z;It((«i,T n u))’ an d (d) average nearest-neighbor 
charge-density correlation 4 \ )), as a function of eaF/tico, for different values of U. The other pa¬ 

rameters are V = 0, g — 0, S = 0.1, a»p = 0, and io = 2.4. In (c), the vicinity of eaF/hcj- 0.9 is enlarged in the 
inset. 


almost independent of U in the range of U < 0.23 [Fig. 4(a)]. The maximum appears around 
eaF/fi(jL)= 0.9. As is clear from the comparison between Figs. 2 and 4(a), it is governed by the 
Bessel function Jo(eaF/ha>) modulating the transfer between sites 1 and 2 within the strongly 
distorted dimer: the system is excited by the electronic transitions within the strongly dis¬ 
torted dimers. The maximum value is indeed close to the energy increment expected when all 
electrons in the strongly distorted dimers transit from the bonding orbitals to the antibonding 
orbitals [Fig. 1(c)]: 2 1\ x 2x(16 sites/4 sites)^ 2.3 x 2 x 4 ^ 18. 

As eaF/ftco increases from zero to about 0.5, the time-averaged charge density ((«,•)) 
changes from the (1,1,2,2) configuration toward a uniform distribution, especially for large 
U, so that the charge order is weakened and the time-averaged double occupancy 
also changes toward a uniform distribution [Fig. 4(b)]. However, around eaF/Tiu>= 0.9, where 
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AZstot reaches the maximum, the configurations of ((«,■)) and almost return to those 

before the photoexcitation, which correspond to eaF/Tico= 0. This indicates that, in the corre¬ 
sponding state, all electrons in the strongly distorted dimers indeed transit from the bonding 
orbitals to the antibonding orbitals [Fig. 1(c)] and almost no electrons are transferred between 
dimers. If only these intradimer transitions were realized, the site-diagonal quantities such as 
((«,)), and the nearest-neighbor charge-density correlation ((m ; m,- + 1 » would be un¬ 

changed. This very state is called a negative-temperature state. 

In the ground state and weakly excited states of small eaFlhio, as U increases, the average 
double occupancy decreases [Fig. 4(c)]. Around eaF/fia)= 0.9, however, the average double 
occupancy for U > 0 is larger than that of the noninteracting ({/=0) ground state (eaF/Twj=0) 
and increases as U increases. This behavior with increasing U corresponds to more attrac¬ 
tive on-site interactions for equilibrium: it appears as if the negative-temperature state has 
an attractive on-site interaction and its strength increases with U when the time-averaged 
correlation is interpreted in terms of the correlation in equilibrium. Near the maximum of 
AE'tot, the energy supplied by the photoexcitation goes slightly to the U term in Eq. (1) to 
enhance the average double occupancy, as discussed later in Sect. 3.5. For small eaF/Tuv, 
as U increases, the average nearest-neighbor charge-density correlation increases [Fig. 4(d)]. 
Around eaF/Tioj=0.9. this correlation is larger than that of the ground state, slightly decreases 
as U varies from U- 0 to (7=0.1, and then increases as U varies from (7=0.1 to U- 0.2. This 
fact cannot be interpreted in terms of an effective on-site interaction for equilibrium. 

The same quantities as in Fig. 4 are shown in Fig. 5 for larger field amplitudes, 3 < 
eaF / Fa) < 8. For comparison, the same ordinate axes are used. Local maxima of the total- 
energy increment AJStot appear around eaF/tuo-4 and eaF/tioj=l3 [Fig. 5(a)]. Their val¬ 
ues depart from the maximum value with increasing eaF/ticj. Around eaF/hcL>= 4, the weak¬ 
ened charge order recovers to a small extent [Fig. 5(b), compare it with Fig. 4(b)]. Around 
eaFI%( jl>= 7.3, the weakened charge order does not recover at all. The realized state differs 
from the typical negative-temperature state at eaF /Ticl>=0.9 because the electrons are trans¬ 
ferred from the weakly distorted dimers to the strongly distorted dimers. Thus, the charge or¬ 
der is weaker and AE tot is smaller than the corresponding quantities at eaF/fia>= 0.9. Around 
eaF/Ti(i)=4 and eaF/fica- 7.3, as U increases, the average double occupancy decreases ex¬ 
cept for eaF/fia>=4 with 0 < U < 0.1 [Fig. 5(c)]: the effective on-site attraction exists for 
eaF4 and 0 < U <0.1 only. The average nearest-neighbor charge-density correlation 
indeed increases as U increases even around them [Fig. 5(d)]. 
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Fig. 5. (Color online) Same quantities as in Fig. 4, with the same ordinate axes and for 3 < eaF/hio < 8. 


3.2 Effect of nearest-neighbor repulsion V 

Time-averaged quantities for U = 0.2 and g - 0 with different values of V and fixed 
co = 2.4 are shown for 0 < eaF/hco < 2.5 in Fig. 6. The total-energy increment A£ lot 
is almost independent of V [Fig. 6(a)]. The maximum appears around eaF/hco= 0.9. The 
maximum value is again close to the energy increment expected when Fig. 1(c) is realized. 
As eaF/hco increases from zero to about 0.5, ((«,}) changes from the (1,1,2,2) configuration 
toward a uniform distribution so that the charge order is weakened and ((n,-n ;+1 )) also changes 
toward a uniform distribution [Fig. 6(b)]. However, around eaF/hco= 0.9, the configurations of 
((«,■)) and ((/?,«,+ 1 )) almost return to those before the photoexcitation. Namely, the negative- 
temperature state [Fig. 1(c)] is realized. 

In the ground state and weakly excited states of small eaF/hco , as V increases, the aver¬ 
age nearest-neighbor charge-density correlation decreases [Fig. 6(c)] and the average double 
occupancy increases [Fig. 6(d)]. Around eaF/hco= 0.9, the average nearest-neighbor charge- 
density correlation is larger than that of the ground state (eaF/hco= 0) for V=0 and increases 
as V increases. At the same place, the average double occupancy decreases as V increases. 





















J. Phys. Soc. Jpn. 




Fig. 6. (Color online) (a) Total-energy increment A£ tot , (b) time-averaged nearest-neighbor charge-density 
correlation ((h,h/ + i)) at interdimer (/ = 2) and two nonequivalent intradimer (/ = 1,3) bonds, (c) average nearest- 
neighbor charge-density correlation 4 YJi=\(( n i n i+ 1))> anc * (d) average double occupancy 4 as a 

function of eaF/Hu), for different values of V. The other parameters are U = 0.2, g - 0, S - 0.1, top - 0, and 
to = 2.4. 


These behaviors with increasing V are consistent with more attractive nearest-neighbor inter¬ 
actions for equilibrium. Near the maximum of AE tot , the supplied energy goes to the V term 
in Eq. (1) to enhance the average nearest-neighbor charge-density correlation, as discussed 
later in Sect. 3.5. 

A typical negative-temperature state does not appear at eaF/fia>=4 or eaF/fico-l .3 (not 
shown), where a local maximum of AEtot appears. Around eaF/T icl>= 4, the weakened charge 
order recovers to a small extent. Around eaF/fito=l.3, the weakened charge order does not 
recover at all. In both cases, as V increases, the average nearest-neighbor charge-density 
correlation decreases except for eaF/Tico=4 with 0 < V < 0.05: the effective nearest-neighbor 
attraction exists for eaF/tiu>=4 and 0 < V < 0.05 only. How these states deviate from the 
typical negative-temperature state is similar to that shown previously in Fig. 5. 
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Fig. 7. (Color online) (a) Total-energy increment AT/ot and (b) time-averaged Holstein electron-phonon cor¬ 
relation averaged over almost singly occupied sites (((( b\ + - 3/2))) + (((£>2 + b\)(ji 2 - 3/2))))/2 and that 

averaged over almost doubly occupied sites ((((£>3 + b\)(iiT, - 3/2))) + (((£>4 + bf)(n 4 - 3/2))»/2, as a function 
of eaF/Hu), for different values of g. The other parameters are U = 0, V = 0. S =0.1, o.>h = 0.2, cjp - 0, and 
w = 2.4. 


3.3 Effect of Holstein electron-phonon coupling g 

The total-energy increment AE tot and the time-averaged Holstein electron-phonon corre¬ 
lation ((( bj + b]){ni - 3/2))) for U = 0, V = 0, and u> H = 0.2 with different values of g and 
fixed oj = 2.4 are shown for 0 < eaF/hco < 2.5 in Figs. 7(a) and 7(b), respectively. The max¬ 
imum of A£ loL appears around eaF/hti>= 0.9 again. The maximum value is close to the energy 
increment expected when Fig. 1(c) is realized: 2fi x 2x(8 sites/4 sites)- 2.3 x 2 x 2 ^ 9. 

In the ground state and weakly excited states of small eaF/ha>, the correlations (((/?, + 
If fn, - 3/2))) are negative and their magnitudes increase as g increases [Fig. 7(b)]. Note 
that, in equilibrium, the magnitude of the correlation ((bj + If )(n, - 3/2)) is larger at the 
almost singly occupied sites even though the intramolecular displacement (bj + if) at the 
almost singly occupied sites and that at the almost doubly occupied sites have opposite signs 
and the same magnitude, as {n t - 3/2) at the almost singly occupied sites (^ -0.5) and 
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that at the almost doubly occupied sites (- +0.5). As eaF/Tico increases from zero to about 
0.7, the correlations are weakened so that the magnitudes of (((£>, + Z?J)(n ( - - 3/2))) become 
small. Around eaF/Tico= 0.9, (((£>,• + b\)(n.j - 3/2))) at the almost singly occupied sites is 
positive and increases as g increases. On the other hand, (((/?, + /+’)(h ( - 3/2))) at the almost 
doubly occupied sites remains negative, but its magnitude is quite small around eaF /fico- 0.9. 
If «(/?,- + b])(n,j - 3/2))) are spatially averaged over all sites, the average correlation is positive 
and increases with g around eaF/boj=0.9. This behavior corresponds to inverted electron- 
phonon couplings (from positive g to negative g) for equilibrium. The energy supplied by the 
photoexcitation goes slightly to the g term in Eq. (1) to invert the average Holstein correlation, 
as discussed later in Sect. 3.5. 

This fact is consistent with the finding in Ref. 15: at the early stage of the photoin- 
duced melting of charge order in a-(BEDT-TTF) 2 l 3 , the oscillation of the correlated electrons 
between the neighboring molecules interferes destructively with intramolecular vibrations. 
Their coupling is theoretically treated by the Holstein term. 15 When a sufficient amount of 
energy goes to the g term, the phase of the electronic oscillation relative to the phase of the 
intramolecular vibration becomes opposite to that in equilibrium. It is because, in the ground 
state, the relative phase is determined in such a way that the interaction energy is lowered. 
This change in the relative phase is observed as the destructive interference. 

A typical negative-temperature state does not appear at eaF/fico= 4 or eaF/tuo-13 (not 
shown), where a local maximum of AEtot appears. Around eaF/hco=4, the magnitudes of 
(((6,+7>J)(«,—3/2))) are small, but they remain negative. Around eaF/%io=l .3, the magnitudes 
of (((/?; + b])(rii - 3/2))) become slightly smaller than the neighborhood, but they are not so 
small compared with those around eaF/Tico=4. Namely, the realized state departs from any 
state with an inverted interaction. 

3.4 Effect ofPeierls electron-lattice coupling S 

Time-averaged quantities for U = 0, V = 0, g = 0, and oo P = 0.05 with different values 
of S and co are shown for 0 < eaF/tioj < 2.5 in Fig. 8. The optical gap sensitively depends 
on S , so that the central frequency co, which is chosen to be resonant with the optical gap, is 
varied: co = 2.19 for S = 0.02, co = 2.25 for S = 0.05, co = 2.41 for S =0.1, and co = 2.8 for 
S = 0.2. As S increases, the energy difference between the bonding and antibonding orbitals 
in the strongly distorted dimer increases. Thus, the maximum of AE tot , which appears around 
eaF/Tico= 0.9, also increases with S [Fig. 8(a)]. The maximum value is again close to the 
energy increment expected when Fig. 1(c) is realized. 
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Fig. 8. (Color online) (a) Total-energy increment A/-’ tot , (b) time-averaged bond density 
((Her ( C I o- C2 >o- + c 2 o- c i,o-)))> ( c ) time-averaged lattice displacement ((u\)), and (d) time-averaged Peierls 
electron-lattice correlation ((u\ Y<r ( c i o- c 2 ,o- + c\ a c \ jCr ))), as a function of eaFItuo, for different values of S. 
The other parameters are U = 0, V — 0, g — 0, and top = 0.05. The central frequency is set at to = 2.19 for 
S = 0.02, to = 2.25 for S = 0.05, to = 2.41 for S = 0.1, and to = 2.8 for S = 0.2. 


As eaF/Tuo increases, the time-averaged bond density between the almost singly occu¬ 
pied sites ((Zcr (c, cr C2,cr + c\ decreases and then becomes negative [Fig. 8(b)]. Around 

eaF/htL>= 0.9, its magnitude shows a local maximum. Although its magnitude is smaller than 
that in the ground state, this corresponds to the negative-temperature state shown in Fig. 1(c). 
Around eaF/htL>= 2, where AE tot almost vanishes, the time-averaged bond density almost 
takes the value taken by the ground state, i.e., the time-averaged bond density almost re¬ 
mains the initial value before the photoexcitation. In the ground state and weakly excited 
states of small eaF/Txto, the time-averaged lattice displacement between the almost singly 
occupied sites ((u\)) is negative and its magnitude increases as S increases [Fig. 8(c)]. As 
eaF/fta) increases, ((wi)) increases and then becomes positive. Around eaF/fict>= 0.9, ((«i)) 
shows a local maximum, which increases with 5. The magnitude of ((i/i)) here is close to, 
but smaller than, that in the ground state. 
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Actually, around eaF/fra>= 2, all time-averaged quantities almost take the values taken 
by the ground state in all cases of Figs. 4, 6, 7, and 8. This situation is similar to the time 
behavior of the double occupancy in the cases eaF/hto=2 and eaF/tioj-l.5 of Ref. 6 after 
the sudden application of a continuous-wave field to the half-filled Hubbard model in the 
dynamical mean-field theory. Although the double occupancy is shown for a short time in 
Ref. 6, its time average would be equal to its initial value for some value of eaF/tiaj slightly 
smaller than 2.5. Note that the dynamical localization is expected for eaF/tica= 2.4 in the 
continuous-wave case. 

In the ground state and weakly excited states of small eaF/tico, the time- 
averaged Peierls electron-lattice correlation between the almost singly occupied sites 
((wi Ho-( c i cr c 2 ,o- + c 2 o- c i>°v» 1S ne g at i ve an d its magnitude increases as S increases 
[Fig. 8(d)]. As eaF/Tica increases from zero to about 0.3, the correlation is weakened so 
that the magnitude of (( U\ Ya (cj ^c^o- + c\ ^Cpo-))) becomes small. Around eaF/hci)= 0.5 and 
eaF/T icl>=\.3, ((w, Y<t ( c ! cr c 2 ,<r + c 2 cr c h<r))) i s quite small but positive for S > 0.05 and in¬ 
creases as S increases. The energy used to invert the Peierls correlation is discussed below in 
Sect. 3.5. The field amplitudes eaF/Tico=0.5 and eaF/fia>= 1.3 correspond to those that weaken 
the charge order to the largest degree in Figs. 4(b) and 6(b). Namely, the electrons are trans¬ 
ferred most extensively between dimers. This is the reason why (( u\ Ya- (cj a c% a + c\ ^c' 1 , 0 -))) 
is quite small. If both the time-averaged bond density and ((u i)) invert their signs, their 
correlation should not invert its sign. This happens around eaF/T ioj- 0.9. Here also, the 
magnitude of their correlation increases as S increases. This corresponds to the negative- 
temperature state. Namely, the negative-temperature state appears when AE tot reaches its 
maximum (eaF/tico= 0.9), while the inverted electron-lattice coupling [(to - uj) to (to + uj) 
in Eq. (1)] appears in different field amplitudes (eaF/1ito= 0.5 and eaF/T ul>=\.3). 

3.5 Energy flows and inverted interactions 

So far, we mainly show whether the correlation functions merely increase or decrease as 
a function of the coupling strength or the field amplitude. In this section, we focus on those 
photoinduced states whose correlation functions behave as if the interactions are inverted. 
We show how the inverted correlations are related to changes in the total energy and in the 
interaction energy. 

The enhanced average double occupancy is defined by 

1 N l N 

— ^ A {/«tt /;T n a » = — ^ (««/,T«U» lt/>o lt/=o) , 


(8) 
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Fig. 9. (Color online) Regions of inverted interactions, denoted by ‘II’, in planes spanned by eaF/tia> and (a) 
U (V=0, g-0, 5=0.1, and cjp- 0), (b) V ((7=0.2, g- 0, 5=0.1, and cj p = 0), (c) g ((/=0, V=0, 5=0.1, u>h= 0.2, 
and cap- 0), and (d) 5 ((/=0, V=0, g=0, and cop- 0.05), from data shown in Figs. 4(c)-5(c), 6(c), 7(b), and 8(d), 
respectively, for resonant excitations. 


for fixed eaF/Tuo. The regions where the enhanced average double occupancy is positive are 
shown as the regions of the inverted interaction in Fig. 9(a). They appear around the local 
maxima of AE tot . For eaF/T ioj=0.92 and £/=0.15, the enhanced average double occupancy 
is 6xl0 -3 in Fig. 4(c). This value is rather small compared with the U dependence of the 
average double occupancy in equilibrium. The energy flowed to enhance the average double 
occupancy AuEp, is defined as 

N 

A u E u = uJ]A u ((n iA n u )) . (9) 

i= 1 

In the case above, A v Eu is 1.4xl0 -2 . We further define A v AE ioX by 

A(jAE lol = AEm lt/>o -AE to t lt/=o , (10) 

for fixed eaFItuo. In the case above, AjjAE tot is 1.8xl0 _1 , while A£ tot is 17. The ratio 
AuEu/AuAEtot is then 7.8x10 2 , while the ratio A f /£ f/ /A£ tot is 8.3x10 '. The increased 
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average double occupancy is defined by 

1 n l N 

— ^ A«n /iT n a » = — ^ («n/,T»a» Ifx> -«»«,T n a)) Ip=o) . (U) 

i— 1 i= 1 

for fixed U. Note that this quantity is negative for U- 0 and positive for U >0.05 in Fig. 4(c). 
Using this quantity, the interaction-energy increment AE U is evaluated as 

N 

AE u = uJ]A((n i ,in ul )). ( 12 ) 

i=i 

For eaF/tico=0.92 and £7=0.15, the increased average double occupancy is 3.2xl0 -2 in 
Fig. 4(c), and A E v is 7.7xl0“ 2 , so that AE V IAE toi is 4.4xl0“ 3 . 

The enhanced average nearest-neighbor charge-density correlation is defined by 

1 N , N 

— ^ Ay«n,n, +1 » = — ^ (««,■«/+1» lv>o -<<«/«,+ 1 » |y=o) , (13) 

i= 1 i=l 

for fixed eaF/Tiio. The regions where the enhanced average nearest-neighbor charge-density 
correlation is positive are shown as the regions of the inverted interaction in Fig. 9(b). They 
appear around the local maxima of AE toX . For eaF/ficj- 0.92 and V=0.2, the enhanced average 
nearest-neighbor charge-density correlation is 3.4xl0 -2 in Fig. 6(c). This value is rather large 
compared with the V dependence of the average nearest-neighbor charge-density correlation 
in equilibrium. The energy flowed to enhance the average nearest-neighbor charge-density 
correlation is evaluated as 

N 

A v E v = vY J ^v((n i n M )) . (14) 

i= 1 

In the case above, A V E V is l.lxlO -1 . The quantity AyAi? to t, defined by 

AyA£ tot = Aiitot |y>o — AF'tot |y=o 5 (15) 

is 5.5xl0 -2 , while A E tot is 17. The ratio A V E V /A v AE tot is then about 2, while the ratio 
A v E v /AE tot is 6.3xl0“ 3 . The fact that A v E v /A v AE tot is larger than 1 is realized by the rela¬ 
tion A v Eu < 0, as shown in Fig. 6(d): the large enhancement of the average nearest-neighbor 
charge-density correlation is realized by the reduction of the average double occupancy in 
the system with competing interactions U and V. The increased average nearest-neighbor 
charge-density correlation, which is defined by 

1 N l N 

— 2 A«n,n (VI » = — 2 («n/n,+i» \f>o -((«,«/+ 1 )) If=o) , 


( 16 ) 
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for fixed V, is 6.5x10 2 in Fig. 6(c), and the interaction-energy increment A E v , 

N 

AE v = vJ]A((n i n i+l )), (17) 

;=i 

is 2.1xl0 _1 , so that AE v /AE tot is 1.2xl0~ 2 in the case above. 

The enhanced average Holstein correlation is given by the average of the quantities for 
the almost singly occupied and doubly occupied sites shown in Fig. 7(b) because it is zero 
for g=0. The region where the enhanced average Holstein correlation is positive is shown as 
the region of the inverted interaction in Fig. 9(c). It appears only around the maximum of 
A E tot . Because we used small g for the technical reason mentioned before, every energy scale 
is small. 

The enhanced average Peierls correlation is given by the quantity shown in Fig. 8(d) 
because it is zero for S =0. The regions where the enhanced average Peierls correlation is 
positive are shown as the regions of the inverted interaction in Fig. 9(d). They appear on both 
sides of the local maxima of A E tot . For eaF/fico=0.51 and S =0.2, this quantity is 1.2x10 2 in 
Fig. 8(d). The energy flowed to enhance the average Peierls correlation A S E S , 

N/2 

AsEs = ^^, ( c 2 /'-1 ,q- C2 .h c r + C 2 /cr^Ty-Po-))) , (18) 

7=1 cr 

is 5xl0“ 2 . The quantity A s A£ tot is defined as before, and A s AE tot is 1.8, while A£ lot is 13. 
The ratio A S E S /A S AE tot is 2.7x10 2 , while the ratio A S E S /AE tot is 3.6xl0“ 3 . The increased 
average Peierls correlation is 7.6x10 1 in Fig. 8(d), and the interaction-energy increment 
A E s , 

N/2 r 

A Es = ^ ((Uj ^ ( C 2j-l,o- C 2j,cr + C 2/,cr C2 7—l.cr))) l^>0 ~(( u j ^ ( C 2;-l,cr C2 7.o' + C 2j,o- C2 j~ l.o - )^ 1^=0 ’ 

j= 1 L (T cr 

(19) 

is about 3.0, so that AE s /AE tot is 0.23. This ratio is relatively large because the ground state 
is stabilized by the Peierls electron-lattice coupling in the first place. 

Comparing Figs. 9(c) with 9(d), we notice that the region of the inverted interaction is 
wider for S than for g. When we see the local maxima of A£ to t- there is indeed a difference 
between them. For g ± 0, with increasing eaF/Tioj. the deviation from the typical negative- 
temperature state becomes large. For S 4- 0 but g = 0, with increasing eaF/fia), the deviation 
does not seem to become so large at least for eaF/Tioj < 10. In the latter, the disordering 
effect seems small. If we increased cl> p and treated Uj quantum-mechanically, the disordering 
effect might become large to suppress the inversion of the Peierls electron-phonon coupling 
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for large eciF/tico. 

In all cases above, for the interaction parameter X (X=U, V, g, or S), the average cor¬ 
relation enhanced by the interaction X is smaller than the average correlation increased by 
the field F. Even in the system with competing interactions U and V, the ratio of the av¬ 
erage nearest-neighbor charge-density correlation enhanced by V to that increased by F 
is about 1 to 2. When the interaction-induced change is compared with the field-induced 
change for the energy ratio, we find A x E x /A x AE tot > AE x /AE tot for X=U and V, while 
A x E x /A x AE l()l < AE x /AE tot for X=S. In any case, A x E x /AE tot is less than 10 -2 and very 
small. The efficiency of the inversion of interactions is quite low. 

4. Discussion 

In general, the formation of a negative-temperature state requires suppressing the dis¬ 
ordering effect (i.e., thermalization) on increasing the total energy. The present model has a 
simple structure. The strongly distorted dimer is regarded as a type of two-level system, which 
is a typical system to achieve a negative-temperature state. The weakly distorted dimers are 
then regarded as the environment, to which the energy flow from the system is allowed. The 
interdimer charge transfer is a type of relaxation process. The structure of the present model 
may be advantageous to forming a negative-temperature state in the sense that the thermaliza¬ 
tion is suppressed by separating the energy scale of the interdimer electron transfer processes 
and that of the intradimer process. 

In real substances, the network of electron transfers is not always simple as treated in 
this study. For instance, in the model previously used for Et 2 Me 2 Sb[Pd(dmit) 2 ] 2, 9 the weakly 
distorted dimer has an unoccupied molecular orbital whose energy is not far from those of 
the antibonding orbitals in the strongly distorted dimer. In such a case, even under optimal 
conditions (such as eaF/ha>= 0.9 in the present study), the electrons can be easily transferred 
between dimers so that the negative-temperature state, if realized, can be short-lived, depend¬ 
ing on the transfer integral between the corresponding molecular orbitals. We expect that, if 
the photoinduced melting of charge order (i.e., a high-temperature state) is easily achieved, 
the negative-temperature state would hardly be achieved, and vice versa. 

5. Conclusions 

Electron-phonon dynamics in the band-insulator phase of one-dimensional three-quarter- 
filled strongly dimerized extended Peierls-Hubbard and Holstein models after the application 
of a symmetric monocycle electric-field pulse with central frequency resonant with the op- 
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tical gap are calculated using the exact diagonalization method. From the field-amplitude 
and interaction-parameter dependences of time-averaged correlation functions, we discuss 
the negative-temperature state and the inversion of interactions. The total-energy increment 
is maximized when the energy gain of an electron transferred to the neighboring site eaF 
under the electric field F is comparable to the photoexcitation energy Tico (eaF/ha>= 0.9). The 
maximum value is close to the energy increment expected when all electrons in the strongly 
distorted dimers transit from the bonding orbitals to the antibonding orbitals, i.e., the en¬ 
ergy difference between the negative-temperature state and the ground state. In this state, 
the site-diagonal quantities such as the time-averaged charge density, double occupancy, and 
nearest-neighbor charge-density correlation are close to those in the ground state, and the 
time-averaged bond density is inverted. Around eaFITicj=4 and eaF/ficj-1 .3, the total-energy 
increment shows a local maximum, but the realized state departs from the typical negative- 
temperature state at eaF/Tuo=0.9. Compared with the latter, the electrons are transferred be¬ 
tween dimers and the charge order is thus weakened. 

The interaction-parameter dependences of the time-averaged correlation functions in the 
negative-temperature state become similar to those in the ground state if the interactions 
(except the Peierls electron-lattice coupling) are inverted: the repulsive on-site (nearest- 
neighbor) interaction is converted to an attractive on-site (nearest-neighbor) interaction, and 
the relative sign between the charge density (n, - 3/2) and the intramolecular displacement 
is inverted. The energy supplied by the photoexcitation goes to these interactions terms. 
As for the Peierls electron-lattice coupling, the relative sign between the bond density and 
the intermolecular displacement is inverted for those field amplitudes (eaF/Tico= 0.5 and 
eaF/tuo=l.3), which are different from the field amplitude for the negative-temperature state 
to appear (eaF/Tico= 0.9). The formation of the negative-temperature state and the inversion 
of interactions have been discussed by Tsuji et al? using the dynamical mean-field theory, 
but they need half-cycle pulses or asymmetric monocycle pulses in the cases studied. In this 
paper, we demonstrate that they are realized even by symmetric monocycle pulses. 
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